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Abstract 

We show that for critical site percolation on the triangular lattice 
two new observables have conformally invariant scaling limits. In par- 
ticular the expected number of clusters separating two pairs of points 
converges to an explicit conformal invariant. Our proof is independent 
of earlier results and SLE techniques, and in principle should provide 
a new approach to establishing conformal invariance of percolation. 

1 Introduction 

Percolation is perhaps the easiest two-dimensional lattice model to formu- 
late, yet it exhibits a very complicated behavior. A number of spectacular 
predictions (unrigorous, but very convincing) appeared in the physics liter- 
ature over the last few decades, see [3] . One of them, the Cardy's formula 
for the scaling limit of crossing probabilities, was recently established for 
the critical site percolation on triangular lattice [IU] . Consequently, scaling 
limits of interfaces were identified with Schramm's SLEq curves, and many 
other predictions were proved, see e.g. |14j . 

In this paper we show that two new observables for the critical site 
percolation on triangular lattice have conformally invariant scaling limits. 
Furthermore, we obtain explicit formulae, consistent with predictions ob- 
tained by physicists [H Our proof is independent of earlier conformal 
invariance results, and uses methods similar to those in [10] rather than 
SLE techniques. It is also restricted to the same triangular lattice model, 
but at least one should be able to use it for a new proof of the conformal 
invariance in this case. 
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2 Notation and Setup 

For convenience reasons, in this paper we shall not work on the triangu- 
lar lattice but rather on its dual, the honeycomb lattice, and thus, rather 
than coloring vertices of triangles, we shall color hexagonal faces (which is 
obviously equivalent). 

2.1 Graph and Model 

Let 11 C C be a Jordan domain (whose boundary is a simple closed curve), 
and orient its boundary dQ counterclockwise. Let I and r be two distinct 
points on dQ, which separate it into a curve u going from r to I (with respect 
to the orientation of d£l) and a curve d going from / to r, such that we have 
dO, \{l,r} = u U d. Let finally w be a point on u. 

Remark 1. The assumption on 0, to be a Jordan domain is not really 
necessary, and the result remains true under weaker assumptions detailed 
in Section 3. We use this assumption in Section^ to avoid lengthy and not 
so interessant discussions. 

We consider the discretization £lg of $7 by regular hexagons defined as 
follows. Let G$ be the regular hexagonal lattice embedded in the complex 
plane with mesh size (i.e. sidelength of a hexagon) 5 > 0. We define Q$ 
as the graph obtained by taking a maximal connected component made 
of hexagonal faces of G$: the union of the closure of the faces is a sim- 
ply connected subset of C. We denote by this subset and by dQ$ the 
(counterclockwise-oriented) simple path consisting in edges of £1$ such that 
£1$ is contained inside it. We define the discretization of /, r and w as the 
closest corresponding vertices of dQ,$, and u and d as the paths from r to I 
and I to r respectively, following the orientation of 80,$. In general we will 
identify I, r, u, d to their respective discretizations. 

We are interested in the process of critical percolation on the faces of 
£1$ : each face of Q$ is colored either in white or in black with probability 
^, independently of the other faces; such a coloring is called a configuration. 
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Figure 1: Domain discretization: the points are marked with crosses, while 
their discretization are marked with circles. The discretization of it and d 
are depicted by sparse and dense dashes respectively. 

More precisely, we are interested in the scaling limit of this model: the 
description of the global geometry of the process as the mesh size 5 tends to 
0. 

Note that for this model \ is known to be the critical value for the 
probability thanks to the work of Wierman and Kesten. However, we do 
not use that this value is critical, only that it is self-dual. 

We call path of hexagons a sequence TL\ , . . . , TL n of hexagons such that 
Til is adjacent to Tii+i f° r i = 1, ■ ■ ■ ,n — 1; a path is called simple if all of 
its hexagons are distinct; a closed simple path (the last hexagon is adjacent 
to the first one) is called a circuit. We say that a hexagon Ti is connected 
to u by a white path if there exists a path of white hexagons that contains 
Ti and that hits u (contains a hexagon having an edge belonging to (the 
discretization of) u). We define similarly connection events involving black 
instead of white paths or connections to d instead of u. We say that a path 
of hexagons 7 separates two families of points A and B if the interior of each 
continuous path a contained in Cl^ from a point of A to a point of B crosses 
the closure of a hexagon of 7. 

We call white cluster a connected (i.e. path-connected in the sense de- 
fined above) set of white hexagons. For a cluster K touching u and d, we 
define its left boundary (respectively right boundary) as the left-most (re- 
spectively right-most) simple path contained in K that touches u and d, i.e. 
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such that there is no path in K separating it from I (respectively r); elemen- 
tary topological considerations show that this notion is well-defined. One 
important property of our lattice is indeed its self- duality: the boundary of 
a white cluster (that does not touch the boundary) is a black circuit, and 
vice versa. 

We will use the term left boundary for simplicity, but strictly speaking 
our definition gives the left-most simple curve inside the cluster, that is 
its left-most boundary after "peninsulas" attached by only one hexagon are 
erased. So this curve would rather bound on the right the dual cluster 
bordering ours on the left. 




Figure 2: In this configuration, the dashed region C is a white cluster, the 
hexagon Ti is connected to u by a white path and the points I and w are 
separated from the points z and r by a black path. 

Notice that since the probability for a hexagon to be white is |, any 
event (i.e. set of configurations) has the same probability as its negative 
with respect to the colors: for instance, the probability that there is a white 
path from x to y is the same as the one that there is a black path from x to 
y. For an event A, we will denote by A the negative event: a configuration oj 
belongs to A if and only if the negative configuration uj (i.e. with the colors 
black and white flipped) belongs to A. 

2.2 Observables 

Let 5 > and consider the process of percolation on as described in the 
previous section. For each z vertex of £1$ we define the following random 
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variables and events: 

• Njj(z): the number of (simple) left boundaries of white clusters touch- 
ing u and d separating I and w from z and r minus the number of 
(simple) left boundaries of white clusters (touching u and d) separat- 
ing I and z from w and r; 

• Ng(z): the same as N l s (z) but for (simple) right boundaries of white 
clusters (also touching u and d); 

• Q'fi(z): the event that there exists a white simple path from d to d 
that separates z from I and r and that is connected to u; 

• Q$(.z): the same event as Qg(z) but with a white simple path path 
from u to u connected to d instead. 

This allows us to define our observables: 



We extend these functions to continuous functions on Clg in the following 
way (in fact any reasonable manner will work): first for the center of a 
hexagon, take the average value of its vertices. Then divide the hexagon 
into six equilateral triangles, and define the functions on each triangle by 
affine interpolation. We can then extend the functions to O in a smooth 
way. 

Remark 2. The point w could in fact be anywhere in £1 (changing its po- 
sition only modifies the functions Hg and Hg by an additive constant). In 
our setup it lies on the boundary for simplicity. 

Remark 3. Another way of writing H l s (similarly for Hg), which motivates 
its definition, is the following: count the expected number of left boundaries 
that separate I from z and r minus the expected number of left boundaries 
that separate I from w and r. 

It is easy to check that this definition is equivalent to the one given above 
(the boundaries that count positively are precisely the ones that separate I 
from r and z but not w, the boundaries that count negatively are the ones 
that separate I from r and w but not z). 

If one uses this way to write H l 5 , taking the difference is essential to get 
a finite limit: as the mesh tends to the expected number of clusters joining 
u to d blows up. 



H\{z) :=E[N l s (z)], 
Hf(z) :=P[QJf(z)], 



Hl{z) :=E[JVJ(z)], 
Hl(z) :=V[Q$(z)]. 
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Remark 4. Notice that the quantities H s and H$ are the same: if one has 
a configuration in {Nl(z) = k}, flipping the colors of all the hexagons gives 
a configuration in {N£(z) = k}. 




Figure 3: In this configuration, 71 and 73 are left boundaries of white 
clusters, 72 and 74 are right boundaries. Q$(z) occurs and we have 
N l s (z) = 0,JVJ(s) = -1 and N\{z') = N r & (z') = 1. 

3 Conformal invariance and main result 

By conformal invariance of critical percolation we mean that the same ob- 
servable on two conformally equivalent Riemann surfaces has the same scal- 
ing limit. 

It was proven in [10] that crossing probabilities of conformal rectangles 
(here the Riemann surface is a simply connected domain with four marked 
boundary points) are conformally invariant and satisfy Cardy's prediction. 

Consequently the interfaces of macroscopic clusters converge to Schramm's 
SLE curves and we can deduce conformal invariance of many other observ- 
ables. 

The goal of this paper is to show conformal invariance of the observables 
+ H£ and H$ — Hf in the same setup, without appealing to the results 
of [TO]. 



6 



3.1 Limit of the observables 

In order to get our conformal invariance result, we prove a more geometrical 
one: a linear combination of our two observables turns out to be (in the 
limit) a conformal mapping. For each S > 0, define H s := {H\ + flj) - 
"§i(Hf - Hf). Then we have: 

Theorem 1. As 5 tends to 0, Hg converges uniformly on the compact subsets 
of£t \ {l, r} to a function h which is the unique conformal mapping from £1 
to the strip S := {x + iy : x € R,y G (— \/3/4, s/3/A)} that maps (in the 
sense of prime ends) I to the left end of the strip , r to the right end of the 
strip and w to ^i. 

Remark 5. The theorem remains valid under the weaker assumption that 
the discretizations £1$ of the domain converge in Caratheodory's sense to 0,, 
in which case the observables converge on the compact subsets ofQ. 

This theorem gives us the asymptotical conformal invariance (and the 
existence of the limit) of the two observables H$ + H$ and H$ — H$ in the 
following sense. 

Corollary 1. Let & be a conformal map as above and denote by H l s , H$, 
Jig and Hf the corresponding observables on the domain Q' := <J?(fi) with 
the corresponding points I' := $(/),r' := &(r),w' := &(w). Then we have 
the following conformal invariance result: 

lim Hl+m = lim (Hi + H r A o $ 

lim Hf - H§ = \im(m - M) o $ 

Proof. By uniqueness of the conformal mapping to S with three points fixed 
we have h = h! o $ (the images of /, r and w by h and h! o $ are the same). 
Taking the real and imaginary parts gives the result. □ 

Taking z and w on the boundary we obtain the conformal invariance of 
the expected number of clusters in a conformal rectangle (a Jordan domain 
with four distinct points on its boundary). Let E be a conformal rectangle 
with the four points 01,02,03,04 in counterclockwise order. Discretize the 
domain and the four points as before and consider the expected number 
C$ of white clusters separating a\ and G14 from 02 and 03, counted in the 
following way: 

• If a cluster touches both (the discretization of) the arcs 0401 and 0203 
(along the counterclockwise orientation of the <9E), it does not count. 
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• If a cluster touches exactly one of the arcs a^ai and 0203, it counts 
once. 

• If it does not touch any of the two arcs, it counts twice. 

Corollary 2. The quantity Cs admits a conformally invariant limit as 5 — > 
0: If S' is another conformal rectangle with the four points a' l5 a' 2 , a 3 , a' 4 , if 
: S — > S' is a conformal mapping such that \P(aj) = a\ for i = 1,2,3,4, 
and is i/ie corresponding number in the domain E', i/ien 

lim C 5 = lim C£ 

Proof. It suffices to take z on the boundary (choose z = ai,w = 02, 1 = 
as, r = 04) and to see that in this case Cs = H l s + Hgi no clusters count neg- 
atively, if a cluster does not touch any arc, both its left and right boundaries 
count, etc. Therefore the result follows from the previous corollary. □ 



3.2 Formulae 

It is not difficult to express the quantity Cs in terms of the cross-ratio (the 
conformal map from the half-plane to a strip is simply a logarithm). If we 
denote by A = ( ai ~ a 3)( a2 ~ Q4 ) the cross-ratio of the four points, we get 

By adding to this formula the probability that a cluster (separating 04 
and a\ from 02 and 03) touches the arc 0401 and the probability that such 
a cluster touches moreover the arc 0203 one can obtain (twice) the expected 
number of clusters without specific counting. 

Using self-duality one can show that these two quantities are the same 
and that they can be expressed as the difference of the probability that 
a cluster separates 04 and a\ from <22 and 03 minus one half times the 
probability that a cluster touches the four sides of our conformal rectangle: 
if there is a black cluster separating 04 and a\ from ai and 03 (event E\ on 
Figure 5), then consider the right-most such cluster; either it touches also 
the arc 0104 (event £3) or it does not; in this latter case by self-duality there 
is a white cluster on its right touching the arcs 0304, 0401 and 0402 (event 
£2). 

Then we can decompose the event £3 in the following way. Either the 
cluster touching the arcs 0304, 0401 and a\Q2 touches also arc 0203 (event 
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£§) or it does not and there is a white cluster that separates it from arc 0203 
(event £4). 

A color-flipping argument gives that events £2 and £4 have the same 
probability (one has that the negative £2 of £2 is £4), which is therefore 
\ (P[fi] - ¥[£&]). Since by self-duality F[£ 3 ] = P[£ 4 ]+P[£ 5 ] we obtain F[£ 3 ] = 
|(P[fi]+P[f s ]). 

Both quantities are conformally invariant and given by Cardy's formula 
(see [3], and [10] for a proof) and by Watts' formula respectively (see [15] . 
and [6] for a proof) . 




Figure 4: Self-duality decomposition 
So one obtains eventually: 
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Proposition 1. The scaling limit of the the expected number of clusters 
separating a\ and 04 from a 2 and 03 is equal to: 
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where the first term comes from Cardy's formula, the second from Watts' 
formula and the third from the main result of our paper. 

3.3 Open questions 

In this paper we show that certain observables have conformally invariant 
scaling limits. The most prominent mathematical tool for rigorous treatment 
of conformal invariance is Schramm's SLE, which describes scaling limits 
of interfaces by the traces of the randomly driven Loewner evolutions - the 
so called SLE K curves, see |8] for an introduction. Once convergence to 
SLE is known, many quantities related to the model can be computed. The 
only proof for percolation uses Cardy's formula for crossing probabilities 
(established for triangular lattice only in |10j ) and locality of percolation or 
the so-called "martingale trick", see [T2"l [S] . 

3.3.1 How to use our observables to establish conformal invari- 
ance of critical percolation? 

Whether our observable can replace the crossing probabilities in the proofs 
above, is interesting even if it has no less dependence on the triangular lat- 
tice. The problem that prevents the direct application of the same technique 
as in [12] is that our observable does not have a "martingale" property (see 
|13j for an overview) with respect to the percolation interface. However, one 
can attempt other approaches, for example exploiting locality. 

3.3.2 Are our observables computable with SLE? 

For the same reason, computing our observables with SLE techniques (using 
this time that the percolation scaling limit is described by SLEq) is not 
immediate. In principle, the computation should be possible, but the setup 
might be difficult. 

3.3.3 Are there other similar observables? 

Similar techniques allow to compute crossing probabilities and two similar 
observables in this paper. One can ask how much more one can learn without 
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appealing to SLE techniques, in particular whether there are any other 
computable observables? 



4 Outline of the proof 

The proof of Theorem [T] consists of three parts. 

• First we prove that from each sequence (H$ n ) n( - N , with 5 n tending to 
0, one can extract a subsequence H$ k which converges uniformly on 
the compact subsets of Q \ {/, r} to a limit function h. 

• We show then that any such subsequential limit h satisfies the following 
boundary conditions: 

Im(h) = on u 

Im(/i) = — on d 
Re(h(w)) = 



• We prove finally that h is analytic. 

In order to get that h is the conformal map (ft of Theorem [TJ we observe 
that h and 4> have the same imaginary part (on the boundary and hence 
inside since the imaginary part is harmonic), and thus have the same real 
part up to a (real) constant by Cauchy-Riemann equations. The constant is 
since the real part of both is equal to at w. Since any subsequential limit 
has the desired value, we conclude by precompactness that Hs converges to 
</>. 



5 Precompactness 

In order to prove the precompactness of the family of functions (Hs)s>o, we 
show that the four families (H l s )$, (H^)s, (Hg)s, (H§)s are uniformly Holder 
continuous on each compact subset of \ {I, r}. Notice that since the inter- 
polation is regular enough we may suppose in the estimates that the points 
we are considering are vertices of the hexagonal faces. 

Lemma 1. For every compact K C £1 \ {/, r} ; the functions and Hf are 
uniformly Holder continuous on K with respect to the metric d of the length 
of the shortest path in fL 
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Proof. We prove the result for Hf. 

Let (3 = m.i z( zK (max (dist(z, u), dist(z, d))). By compactness of K we 
have (3 > 0, and so each point in K is at distance at least (3 from u or d. 

We have that for each z, z' G r2 the disc D := L>((z + z')/2,d(z, z')) 
contains a path from z to z' . 

Since \H^\ is uniformly bounded, we can assume from now that the points 
z and z' (in K) are close enough, i.e. such that d(z, z') < (3. By elementary 
partitioning we have that \H$(z) - H%(z')\ < P[QJf(z) \ Qs( z ')} +^[Qs(z') \ 
Q'^z)]. So it is enough to show that there exists C > and a > such that 

nQ u 5 {z)\Q u s {z')]<C .d(z,zY. 




Figure 5: The event Q^{z) \ Q^{z') implies a black connection of a micro- 
scopic circle of radius £(\z — z'\) to a macroscopic circle of radius (K is 
not fixed precisely on this picture). 

By self-duality, we have that the occurence of the event Q^{z) \ Q^(z') 
implies the connection of the boundary of the disc D to u by a black path 
and to d by two disjoint white paths. Since at least one of the two sides is 
at distance (3 (for 5 sufficiently small, which we may suppose), this event 
implies the connection (by a black or white path) of a (microscopic) circle 
of radius d(z,z') to a circle of (macroscopic) radius (3. By Russo-Seymour- 
Welsh Theorem (see [2], [7] for instance), there exists C > and a > such 
that this event is of probability less that C ■ d(z, z') a (uniformly in 5) and 
this gives us the desired result. □ 
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Lemma 2. For every compact K C O \ {I, r}, the functions H s and HI are 
uniformly bounded and uniformly Holder continuous on K with respect the 
metric d of the previous lemma. 

Proof. The proof is essentially the same as for the previous lemma: the 
probability that a cluster passes between two close points z and z' is small 
(say C(z, z')) for the same reasons. To control the expectation, we can use 
BK inequality which gives that the probability that n disjoint clusters pass 
between z and z' is smaller that C(z, z') n . 



Figure 6: The event {Ng(z) = Nl(z') + 1} implies a white connection of a 
microscopic circle of diameter d(z, z') to a macroscopic circle of diameter (3 
(K is not fixed precisely on this picture). 



Proposition 2. The function family (Hs)s>o is precompact with respect to 
the topology of uniform convergence on every compact subset of£l \ {l,r}. 

Proof. We are only interested in letting 5 tend to (and otherwise it is any- 
way trivial). So let 5 n be a sequence tending to 0. On each compact subset 
if of O \ {l,r}, the functions H l s ,H^,H^,H^ are bounded and uniformly 
Holder continuous in 5, so they form equicontinuous families. By Arzela- 
Ascoli's theorem, they form a precompact family. We can therefore extract 
a subsequence 5^ of 5 n such that H\ k , Hg , Hf , H$ converge uniformly on 
K. Since Cl \ {I, r} can be written as a countable union of compact subsets, 
a diagonal extraction gives us the desired result. □ 




U - 



□ 
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6 Boundary conditions 



Lemma 3. We have the following boundary conditions: 



8->o 



lim H%(z) = 0, 




Mz G u 




lim Hfiz) = 0, 



lim mtw) = 



Vzed 



Proof. By definition and continuity the condition for H s and Hg is obvious. 

For the first boundary value, notice that for z on u, the event Qg(z) 
implies the connection of z to d (which is at a positive distance from z) by 
two white paths. By Russo- Seymour- Welsh, this probability tends to as 
5 — > 0, so we are done. 

For the second one, first notice that both Qg(z) and its color-negative 

Qx(z) cannot occur simultaneously for z G u. By symmetry we obtain that 
H${z) < K- To see that the limit is actually I, it suffices because of the 

symmetry to prove that the probability that neither Qg(z) nor Qg(z) occur 
tends to as 5 — > 0. 

Indeed if Q$(z) does not occur, then either there is no black path separat- 
ing z from d (call this event A) or there is at least one black path separating 
z from d but these black paths do not touch d (event B). By self-duality, 
A is the event that z is connected to d by a white path. Again by self- 
duality, the occurence of B implies that Qi{z) occurs: take the lowest black 
path 7 separating z from d (which does not touch d by definition), so its 
lower boundary is a white path that touches d (otherwise this white path 
would have a lower boundary which would be a black path and would thus 
contradict the definition of 7), which implies that Qg(z) occurs. 

So if neither Qg(z) nor Qg(z) happen, A happens. But as seen above, 
the probability of A tends to 0, since the probability of a connection by a 
white path from z to d tends to as 5 — > 0. 

The arguments for z £ d are the same as the ones for z G u. 



7 Analyticity 

We are now interested in showing the analyticity of any subsequential limit of 
Hg = (H l g + H r s ) - ^§-i(Hf - Hf) as 5 -» (since by Proposition the family 



□ 
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Figure 7: With probability tending to one as 5 — > 0, exactly one of these 
two events occurs. 

of functions (Hg) s>0 is precompact). The main step consists in proving that 
for each 5 > 0, the function Hg is discrete analytic in a sense explained in the 
next paragraph, which allows to show that Morera's condition is satisfied. 

7.1 Discrete Cauchy-Riemann equations 

Let us first introduce several notations. 

For an oriented edge e = {x, y) in the interior of Qg, let us denote by rTe 
and T 2 .e the edges of Qg obtained by rotating counterclockwise e around x 
by an angle of 2tt/3 and 4tt/3 respectively. We will denote by e* the dual 
edge of e: the edge from the center of the hexagon on the right of e to the 
center of the hexagon on the left of e. 

For a function F defined on the set of vertices of ttg and an oriented edge 
e = (x,y), let us define dgF as F(y) — F(x). 

Let dgH\ be P[N l s {y) = N l 5 (x) ± 1] By linearity of the expectation it 
is easy to see that d s H\ = d±H\ - dZ H\. Let d$H% be F[C}f(y) \ Q u 5 {x)} 
and dZHI be d±gHf = P[QJf(x) \ QJf(y)]. As before we have dgHf = 

For H r s and Hg, we define dJ and dg in the same way as for H l s and 
respectively and also obtain dg = di — d~ . By linearity it is also defined for 
Hg. 

We have the following discrete analyticity result, which already suggests 
that Hg is analytic in the limit and is a discrete analogue of the Cauchy- 
Riemann equations. 

Proposition 3 (Discrete Cauchy-Riemann equations). For any S > and 
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Figure 8: Edges notations 



any oriented edge e in the interior oftts, we have the following identity: 




Proof. Notice that since each configuration (coloring of the hexagons) has 
equal probability, bijective maps are measure-preserving. We will use this 
fact several times in the proof. Fix 5 > 0, take as before e = {x, y) and 
introduce the following notations. 

In what follows, r.y and r 2 .y will be the vertices of Qg such that rTe = 
(x,ry) and r 2 .e = (x,T 2 y). Let C (respectively 1Z ; Z ; T) be the hexagonal 
face that is adjacent to e and r 2 .e (respectively to e and r7e ; to rTe and 
r 2 .e ; the hexagon that touches y). For a hexagonal face, for instance C, we 
denote by C w the event that this face is connected by a white path to d, by 
C b the event that it is connected by a black path to u, by the event that 
it is connected by (not necessarily disjoint) white paths to both u and d, 
and etc.: the connections to u are denoted by superscripts, the connections 
to d by subscripts. Recall that we use the notation Ao B for the event that 
both A and B occur on disjoint sites (notice that it is well defined for the 
events we use here). 

We now compute the derivative dZ of Hg. We have that the event 
A := {N%(x) = N T 5 (y) + 1} is the same as B := X\ o C w o R w , since it is 
clear that B implies A, and by self-duality, if B does not occur, then A does 
not occur (since otherwise there would be a white path touching the right 
boundary of the white cluster passing between y and x and separating x 
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from r which would be absurd by definition of the right boundary) , so both 
are equal. 

Notice that on this event, by going from y to x, either we gain a clus- 
ter boundary counting positively or we lose a cluster boundary counting 
negatively 

If B occurs, then we can define A as the counterclockwise-most extremal 
white path that joins C to d (call A^ its hexagon on d) and p as the clockwise- 
most extremal white path that joins 1Z to u (call p u its hexagon on u). We 
can then us a self-duality argument in the interior of rectangle I, Xd, T, p u (we 
consider the topological rectangle delimited by A (excluded), p (excluded), 
the arc p u l (included) and the arc (included)): B is the disjoint union 
of C and D, where C is the event that B happens and that there is a white 
path that joins the arcs l\d and Tp u and D is the event that B happens 
and that there is a black path that joins the arcs p u l and A^T (these events 
occur in the interior of the rectangle). So we have ¥[B] = P[C] + F[D]. But 
C is equal to X\ o C w o and we have that D and X\ o o TV" are clearly 
in bijection: it suffices to flip (i.e. invert) the colors inside the rectangle to 
map one onto the other (this is well-defined because the definition of the 
rectangle does not depend on the colors of the hexagons inside) , and so the 
configuration inside is independent of the colors elsewhere. 

But now we have that X b o C w o TV£ and X™o£ b o TV£ also have the same 
probability. Let i be the clockwise-most extremal black path that joins X to 
u, and flip the colors in the interior of the part of the graph G comprised 
between i and A that contains I (i and A excluded). Then flip all the colors of 
£ls- This defines a (clearly bijective) map from X\ oC w oTZZ to X™o£ b olZ™. 
The same color-flipping argument shows that X\ o C™ o TV" and X™ o £™ o TZ b 
also have the same probability. So we can summarize the discussion above 
in the following equations, see Figure 10: 

dZHl = ¥[Nl{x) = Nl{y) + l] 

= F[X b b oC w oTV"} 

= ¥[X b b o C w o TV" W \ + P[J 5 6 o CI o TV"} 

= P[X>£ 6 o7C]+P[J>£>^]. 

Using a very similar method (but considering this time a rectangle that 
contains r instead of I when applying self-duality), one obtains, see Figure 
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11: 



d+H l s = F[N l s (y) = N l s (x) + l] 
= W[l h b o£ w oK w } 

= ¥[i b b ocZon w \+¥[i b b oc b onz\ 




Figure 9: Computation of d- H$. White paths are dashed and black path 
bold. The stripped region is G (the subgraph where the color are flipped). 
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Figure 10: Computation of dtH\. 



Let us now compute the derivative c^t g of Hg. By self-duality we have 
that the event X := Qg(r.y) \ Qg(x) is the same as the event that I and 
C are on a white simple path from d to d which is connected to u. 1Z is 
connected by a black path to u (otherwise there would be a white path 
separating it from u and this path would be connected by a white path to u 
as well because Q^(r.y) occurs, which would imply that Q^{x) also occurs). 
Suppose that X occurs. Let A' be the clockwise-most extremal white path 
that joins C to d and i! the counterclockwise-most extremal white path that 
joins X to d. Then obviously, exactly one of the three following events occurs: 



19 



1. Y := T™ o £™ o 1Z b : there is a white path that joins A' to u and there 
is a white path that joins i! to u. 

2. Z :=Z™ o C w o 1Z b b : there is a white path that joins A' to u and there 
is no white path that joins i! to u. 

3. W := 1 W o C™ o 7?|: there is no white path that joins A' to u and there 
is a white path that joins i' to u. 

Using a color- nipping argument we obtain that F[Z] = P[Z^ o Cb ° 7?.™]: 
take the counterclockwise-most black path that joins 1Z to u, call it 71, the 
clockwise-most white path that joins C to d, call it 72, flip the colors in the 
interior of the part of £1$ delimited by 71 and 72 that contains r (71 and 72 
excluded), then flip the colors of the whole graph. This defines a bijection 
from Z to 1™ o Cb o 72.™ . Thus we have shown that 

a+ e H2 = pko£:o^ip[i:oaok:]ip[v£:ok:]. 

One obtains similarly, see Figures 13-15 at the end of the section: 

d+H$ = P[J>£ 6 o7C]+P[J>£>7e 6 ]+P[J 6 o£ 
8+ H% = PR C K b ] + PR o £ 6 o ft™] + P[J 6 o £ 

T .e 

9+ i# = PR o L h o ft™] + PR C + ^ 6 £ 

r .e 

Summing up the identities obtained so far, we obtain the desired result. □ 
7.2 Morera's condition 

The last step in order to prove the analyticity of h is to show that any 
contour integral of the subsequential limit h vanishes. This is given by the 
following proposition (since the convergence is uniform on each compact 
subset of S7, the integral is equal to the limit of the integrals § H$(z) dz as 
5^0). 

Proposition 4 (Morera's condition for h). Let 7 be a simple closed smooth 
curve in Q, oriented counterclockwise. Then we have 

H s (z) dz^O as 5 -> 

Proof. For each sufficiently small 5 > 0, let 75 be a discretization of 7, such 
that 7,5 is a simple curve oriented in the same direction consisting in edges 
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that follow the orientation of y$ and such that 75 —* 7 as (5 — > in the 
Hausdorff metric and with a number of edges of order 

For e = (x,y), let us define F(e) := ^M+jw anc i g* = y _ x ( w hen 
appears alone). We approximate the integral <f^Hs(z)dz by a Riemann 
sum along 7,5 defined as ^2g &y ^s(^)- 



As <5 — > 0, one has 



0, by precompactness 



f 7 ^)d^-E ?e7 e^(e) 
of the family (Hg)g >0 in the topology of uniform convergence on the compact 
subsets. 

We now use the following discrete summation lemma (cf. [1]). Define 
Int(7,$) as the set of all oriented edges lying in the interior of the part of 
which is inside 75 and recall that e* is the dual edge of e (seen as a scalar it 
is equal to y/3ie). 

Lemma 4. 

Y j eH 8 {e)= e*d+H s +o 5 ^ (l) 

e675 eelnt(7«) 

Proof. Denote by Hex^) the set of hexagonal faces of £1$ which are inside 
75 and for such a face /, denote by df the set of its six edges oriented in 
counterclockwise direction. We have that 

ee 7 * /eHex( 75 )/eO/ 

since the terms appearing in edges that are not on 75 appear twice (in 
two faces to which such an edge belongs) with opposite signs and therefore 
cancel. Denote by {xo, xi), (x\, X2), ■ ■ ■ , (#5, Xq) the six edges of df and take 
the indices modulo 6; denote by c(/) the center of a hexagonal face (this 
term is purely artificial yet). A simple calculation shows: 

£ eH s (e) = £ ( Xfc+ / fc+1 - c(f)) (H(x k+l ) - H(x k )). 

eed/ fc=0 ^ ' 

If {xk,x k + 1) does not lie on 75, the term Xk+Xk + 1 (g(x fc+1 ) — H(xk)) 
appears twice with opposite signs and cancels, so only the terms with the 
factor c(/) remain. A term of the form H{x k+ i) — H{xk) becomes a factor of 
the difference between two center faces which is the edge dual to (x^, Xk+i}- 

On the other hand, we have that the contribution of the boundary terms 
on 7,5 tends to 0: we have that the number of edges of 75 is of order 5~ 1 , the 
term Xk+Xk + 1 _ C (J) j s Q f order 5 and H$ is Holder on a neighborhood of 7. 
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We obtain that the sum is equal to 

e*d g H s + o(l), as 5^0 

eeCcwInt(7(s) 

where Ccwlnt is the set of the counterclockwise oriented edges of the set 
of faces Hex(7,5). Taking the sum over the set Int^a) of all oriented edges 
inside 7,5, using dg = dt — <9+-, we obtain 

^2 e*dtH s + o(l), as 5^0 

eelnt(7 5 ) 

as required. 

□ 

Now it suffices to prove that the sum ^eeint^) e *dgHg given by the 
previous lemma is equal to 0. This is given by the discrete Cauchy-Riemann 
equations. Let us reorder the terms in the sum in the following way: 

£ e*d+H s = £ e*d+ ( H r s + H l s - ^ i(Hf - H$) ) 

e6lnt(7 5 ) eelnt^) V / 



- Y, e*(dZHl-dtH\ + ^idt{H%-H$)\ 

e6lnt(7^) V / 



'd+m - a+H$ - d+ m + 9- h 



eelnt(7j) 

= - V e* 

eelnt(7 5 ) \ 

where last equality is obtained using the discrete Cauchy-Riemann equa- 
tions. Reordering one last time the sum (using the changes of variables 
(r.e)* — ► e* and (r 2 .e)* — > e* in the first and the second parts of the sum 
respectively), we obtain 



2 

eelnt(7 4 ) 



which is equal to 0, since yy?>ie* + ((r 2 .e)* — (r.e)*)J = by the geometry 

of the lattice (and this is in fact the only step in our proof where the actual 
embedding of the lattice is crucial) . 

□ 
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Figure 11: Computation of d'uHg. 
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Figure 12: Computation of Hf . 
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Figure 13: Computation of H%. 

-7-2 p 
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Figure 14: Computation of d + ~. H$. 
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